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Abstract 

This  paper  shows  the  importance  of  including  the  activation  losses  in  a  model  designed  to  accurately  predict  the  voltage-current  relationship 
of  a  modem  solid  oxide  fuel  cell.  The  governing  equation  for  voltage  drop  due  to  activation  losses,  the  Butler- Volmer  equation,  is  discussed 
and  the  derivation  and  applicability  of  several  approximations  is  investigated.  Finally,  the  assumptions  required  to  combine  the  separate  anode 
and  cathode  activation  loss  calculations  into  a  single-term  equation  are  clarified. 
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1.  Introduction 

The  potential  promise  of  electricity  generation  with  higher 
efficiency,  lower  cost,  and  lower  pollution  using  solid  oxide 
fuel  cell  (SOFC)  power  systems  has  fuelled  an  unprecedented 
wave  of  research  focused  on  developing  and  commercializ¬ 
ing  this  technology.  Computational  models  are  being  devel¬ 
oped  by  many  research  groups  [1-17]  to  study  the  behavior 
of  complex  SOFC  power  systems  and  to  help  direct  efforts 
for  increasing  stack  performance  and  lowering  installation 
costs.  The  level  of  detail  and  applicability  of  these  models 
depends  on  the  intended  purpose.  Some  are  used  for  high 
level  design-point  and  off-design  system  scoping  studies, 
while  others  are  focused  on  system  behavior  during  tran¬ 
sient  operation.  Lower  level  models  are  developed  to  better 
understand  detailed  reaction  mechanisms  within  each  cell  to 
guide  improvements  in  materials  selection  and  stack  manu¬ 
facturing  processes.  Regardless  of  the  intended  application, 
a  common  goal  is  necessary  to  accurately  characterize  the 
voltage-current  relationship  of  the  SOFC  unit. 
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When  the  SOFC  is  not  connected  to  an  external  load,  there 
is  no  current  flow  outside  of  the  stack  and  the  operating  volt¬ 
age  is  equal  to  the  open-circuit  voltage.  The  open-circuit  volt¬ 
age  for  an  individual  cell  is  typically  assumed  to  be  equivalent 
to  the  equilibrium  Nernst  potential  for  the  overall  reaction 
H2  +  l/202  -*  H20. 


This  assumes  any  irreversible  processes  occurring  at  the 
open-circuit  condition  are  negligible,  such  as  reactant 
and  electronic  electrolyte  crossover  (i.e.,  internal  short- 
circuiting)  [1,18]  and  electrode-electrolyte  parasitic  reac¬ 
tions  (e.g.,  corrosion-type  reactions  that  can  significantly  re¬ 
duce  the  open-circuit  voltage  of  other  types  of  fuel  cells) 
[2,19]. 

As  the  stack  begins  to  supply  current  to  an  external 
load,  it  moves  away  from  a  state  of  equilibrium  and  the 
operating  voltage  drops  due  to  irreversibilities  associated 
with  internal  charge  transfer,  conduction,  and  diffusion 
processes.  These  irreversibilities  are  separately  categorized 
as  activation,  ohmic,  and  concentration  losses  (it  should 
be  noted  the  term  “loss”  used  here  is  equivalent  to  the 
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Nomenclature 

A 

pre-exponential  used  in  the  exchange  current 
density  equation  [A  cm-2] 

£act 

activation  energy  [J  mol-1] 

F 

Faraday’s  constant  [96,485  C  mol-1  or 
Jmol-1  V-1] 

AG 

change  in  Gibbs  free  energy  [J  mol-1] 

A  G° 

change  in  Gibbs  free  energy  at  1  bar  [J  mol-1] 

i 

current  density  [A  cm-2] 

z'o 

exchange  current  density  [A  cm-2] 

ZD, A 

anode  exchange  current  density  [A  cm-2] 

*'o,c 

cathode  exchange  current  density  [A  cm-2] 

ZO.eff 

effective  exchange  current  density  [A cm-2] 

z'l 

limiting  current  density  [A  cm-2] 

ne 

number  of  electrons  transferred  in  reaction 

Pi 

partial  pressure  of  gas  constituent  z  [bar] 

Ru 

universal  gas  constant  [8.3 14  J  mol-1  K-1] 

T 

temperature  [K] 

^Nemst 

equilibrium  Nernst  potential  [V] 

Voc 

open-circuit  voltage  [V] 

Voperatin; 

g  cell  operating  voltage  [V] 

AVact 

activation  loss  [V] 

AVact,  A 

anode  activation  loss  [V] 

AVact'c 

cathode  activation  loss  [V] 

A  Vact,totai  total  cell  activation  loss  [V] 

AVconc 

concentration  loss  [V] 

AVohmic 

:  ohmic  loss  [V] 

Greek  letters 

Oil 

reduction  charge  transfer  coefficient 

at, A 

anode  reduction  charge  transfer  coefficient 

at,c 

cathode  reduction  charge  transfer  coefficient 

Oil 

oxidation  charge  transfer  coefficient 

«eff 

effective  charge  transfer  coefficient 

P 

charge  transfer  symmetry  factor 

terms  “overpotential”  and  “polarization”  commonly  found 
in  electrochemistry  literature).  Activation  losses  are  as¬ 
sociated  with  overcoming  reaction  energy  barriers  at  the 
electrode-electrolyte  interfaces.  Ohmic  losses  are  associated 
with  electron  and  ion  conduction  processes  occurring  in 
the  electrodes,  electrolyte,  and  interconnects,  as  well  as 
the  contact  resistances  across  each  material  interface.  Con¬ 
centration  losses  are  associated  with  reactant  and  product 
diffusion  limitations  between  the  bulk  flow  and  reaction 
sites. 

^operating  =  V0C  -  AVact  -  AVohmic  -  AVconc  (2) 

The  purpose  of  this  paper  is  to  compare  several  meth¬ 
ods  found  in  literature  for  calculating  the  activation  loss 
and  to  clarify  some  common  misconceptions.  Methods 
for  calculating  the  ohmic  and  concentration  losses  are 
equally  important  and  highly  dependent  on  cell  geome¬ 


try  and  manufacturing  processes.  The  reader  is  directed 
to  several  good  references  on  these  subjects  [1-17, 
20-23], 

A  common  assumption  made  in  the  past  SOFC  models  was 
to  neglect  the  activation  loss  all  together.  Two  main  reasons 
typically  given  were:  (a)  SOFCs  operating  around  1000  °C 
have  very  rapid  chemical  kinetics,  which  tends  to  signifi¬ 
cantly  lower  activation  losses,  and  (b)  the  ohmic  losses  are 
much  greater  in  magnitude  and  dominate  the  overall  cell 
behavior  under  normal  operating  current  densities.  These 
assumptions  may  be  valid  for  the  older  tubular  geometries 
which  were  originally  designed  to  operate  at  this  high  tem¬ 
perature  and  had  relatively  long  current  paths,  but  modem 
stacks  are  being  designed  to  operate  at,  or  below,  800  °C 
and  utilize  planar  geometries  with  very  short  conduction 
paths.  This  results  in  activation  losses  closer  in  magnitude 
to,  or  even  greater  than,  the  other  two  loss  mechanisms  (see 
Fig.  1). 

Another  potential  pitfall  for  modelers  is  the  common  mis¬ 
conception  that  the  predominantly  linear  profile  of  many 
SOFC  voltage-current  curves  is  the  result  of  ohmic  loss  dom¬ 
ination.  Since  ohmic  losses  are  a  linear  function  of  current 
density  and  both  activation  and  concentration  losses  are  log¬ 
arithmic  functions  of  current  density,  the  danger  is  to  assume 
these  latter  two  losses  do  not  significantly  contribute  to  the 
overall  voltage  drop  when  the  voltage-current  curve  is  nearly 
linear. 


AVohmic  0^  i 


(3) 


(4) 


AVconc  OC  In  1  - 


(5) 


where  z'o  (equilibrium  exchange  current  density)  is  defined 
as  the  equal  forward  and  reverse  flow  of  electrons  across 
the  electrode-electrolyte  interface  under  open-circuit  con¬ 
ditions  and  /l  (limiting  current  density)  is  defined  as  the 
maximum  possible  current  density,  in  which  the  reactants 
are  consumed  at  their  maximum  possible  replenishment 
rates. 

In  reality,  the  activation  losses  tend  to  steadily  increase 
over  the  operating  range  of  the  SOFC  due  to  relatively  large 
values  of  the  equilibrium  exchange  current  density,  z'o.  In 
contrast,  low  temperature  fuel  cells  tend  to  have  smaller 
exchange  current  densities  and  consequently  a  rapid  in¬ 
crease  in  activation  losses  at  low  operating  currents  fol¬ 
lowed  by  a  leveling  out  at  higher  currents.  Also  evident 
in  Fig.  1  is  the  relatively  small  magnitude  of  the  con¬ 
centration  loss  over  most  of  the  operating  range,  thus  it 
does  not  significantly  affect  the  slope  of  the  operating 
voltage  curve  until  the  limiting  current  density,  z'l,  is  ap¬ 
proached. 
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2.  Butler- Volmer  equation  and  approximations 

The  governing  equation  for  the  activation  loss  at  each 
electrode-electrolyte  interface  (i.e.,  the  equation  is  ap¬ 
plied  separately  at  the  anode  and  cathode)  is  the  general 
Butler- Volmer  (B-V)  equation  [19,22]: 
i  =  i0(e“Xf/SuT) av«  _  e-«2(F/wy act)  (6) 

The  equilibrium  exchange  current  density,  to,  is  a  strong  func¬ 
tion  of  cell  materials,  construction,  and  temperature  and  is 
typically  calculated  using  an  Arrhenius-type  relation  of  the 
form: 

i0  =  Ae-E^R“T  (7) 

where  the  pre-exponential  term.  A,  is  generally  a  function  of 
cell  temperature  and  reactant  and  product  partial  pressures 
[1-4, 6, 8, 9,1 1,15].  The  reduction  and  oxidation  transfer  coef¬ 
ficients,  ai  and  a2  respectively,  are  governed  by  the  electron 
transfer  processes  occurring  across  the  electrode-electrolyte 
interface.  The  net  reactions  at  the  anode  and  cathode  are 
well  understood  for  a  hydrogen  fueled  SOFC  and  given 
by  [H2  +  O2-  -*  H20  +  2e]  and  [l/202  +  2e  -»  02“],  respec¬ 
tively.  There  may  be  several  intermediate  steps,  however,  and 
the  exact  reaction  mechanisms  are  difficult  to  identify.  If  the 
mechanisms  are  known,  the  published  method  [19,22]  for 
determining  the  transfer  coefficients  is  as  follows: 


occurs  for  one  act  of  the  overall  reaction,  s  the  number  of 
electrons  transferred  in  the  rds  (0  or  1)  and  f)  the  symmetry 
factor  (defined  below). 

A  simplification  made  in  many  SOFC  models  is  to  assume 
each  reaction  occurrence  is  a  one-step,  single-electron  trans¬ 
fer  process.  This  results  in  the  following  form  of  the  B-V 
equation: 

I  =  ^(eO-AK^uOAV*,  _  g—f(F/RuT)A vact •,  (10) 

where  the  symmetry  factor,  fi,  is  defined  as  the  fraction  of 
the  activation  voltage  loss  that  affects  the  activation  energy 
barrier,  and  thus  the  rate  of  electrochemical  transformation. 
For  simplicity,  this  factor  is  commonly  assumed  to  be  0.5  for 
fuel  cells,  but  values  ranging  from  0.3  to  0.6  have  been  used 
to  match  experimental  data  [9]. 

Since  the  full  B-V  equation  must  be  solved  implicitly  for 
the  activation  voltage  loss,  several  explicit  approximations 
have  been  utilized  in  the  literature.  It  is  crucial  to  under¬ 
stand  the  range  of  applicability  of  each  of  these  equations 
to  minimize  modeling  inaccuracies.When  the  activation  loss 
is  large  (AVact>200mV),  the  first  exponential  term  in  the 
B-V  equation  is  much  larger  in  magnitude  than  the  second. 
If  the  second  exponential  is  neglected,  the  resulting  equation 
is  the  high-field  approximation  (or  sometimes  called  the  Tafel 
equation  [1,7,13,18,19,21-23]).  This  equation  is  valid  when 
j/t’o  >  4,  as  discussed  in  the  next  section: 


«i  =  -  sfi  (8) 

a2  =  -  +  sfi  (9) 

where  n  is  the  total  number  of  electrons  transferred  in  the 
overall  reaction,  y  the  number  of  electrons  transferred  before 
the  rate  determining  step  (rds),  v  the  number  of  times  the  rds 


oc\F  \ioJ 


(11) 


If  the  exponential  terms  in  the  one-step,  single-electron  trans¬ 
fer  process  B-V  equation  are  expanded  in  a  power  series  and 
the  higher-order,  non-linear  terms  are  neglected,  the  resulting 
equation  is  the  low-field  approximation  (or  sometimes  called 
the  linear  current— potential  equation  [1,7,13,19,21-23]). 
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This  equation  is  valid  when  i/io  <  1,  as  discussed  in  the  next 
section: 


Another  approximation  to  the  B-V  equation  comes  from  rec¬ 
ognizing  the  shape  of  the  corresponding  A  Tact  versus  current 
curve  is  similar  to  a  hyperbolic  sine  function.  In  fact  if  a\  and 
0J2  are  assumed  to  be  equal,  the  activation  loss  can  be  explic¬ 
itly  found  using  [19,22]: 


This  will  be  referred  to  as  the  hyperbolic  sine  approximation. 


3.  Applicability  of  approximations 

As  stated  earlier,  the  applicability  of  each  explicit  approx¬ 
imation  to  the  implicit  B-V  equation  must  be  analyzed  to 
avoid  significant  modeling  errors.  To  simplify  the  analyses,  a 
one-step,  single-electron  transfer  process  is  assumed  and  the 
current  density  is  normalized  by  the  equilibrium  exchange 
current  density.  Fig.  2  shows  a  comparison  of  activation  loss 
calculations  for  a  cell  operating  at  800  °C  and  a  symmetry 
factor  of  0.5.  Note  the  hyperbolic  sine  approximation  is 
exactly  equal  to  the  B-V  equation  for  this  case.  The  low-field 
approximation  is  only  applicable  when  the  operating  current 
density  is  less  than  or  equal  to  the  equilibrium  exchange 
current  density.  The  high-field  approximation  is  only 


i/i0 

Fig.  2.  Comparison  of  high-field  and  low-field  approximations  to  the  Butler-Volmer  equation  for  a  one-step,  single-electron  transfer  process. 


current  density  (A/cm2) 


Fig.  3.  Comparison  of  experimental  voltage-current  data  [9]  and  model-t 


.  curves  using  different  activation  loss  equations. 
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applicable  when  the  operating  current  density  is  greater 
than  four  times  the  equilibrium  exchange  current  density. 
There  is  a  significant  gap  (l<i/io<4)  where  neither  of 
these  approximations  adequately  determines  the  activation 
loss. 

Fig.  3  shows  the  applicability  of  each  activation  loss 
equation  used  in  a  cell  model  validated  with  experimen¬ 
tal  data  from  a  modem  planar  SOFC.  The  model  results 
using  the  B-V  equation  and  hyperbolic  sine  approxima¬ 
tion  are  nearly  identical  and  match  well  with  the  exper¬ 
imental  data.  The  model  results  using  the  high-field  and 
low-held  approximations  deviate  significantly  from  the  ex¬ 
perimental  data,  except  in  very  limited  ranges  of  current 
density. 


To  better  understand  the  validity  of  each  approximation 
relative  to  the  B-V  equation  as  a  function  of  symmetry  fac¬ 
tor,  an  analysis  was  performed  where  the  activation  loss  er¬ 
ror  is  calculated  over  a  wide  range  of  current  density  ra¬ 
tios  and  symmetry  factors.  Fig.  4  illustrates  the  error  us¬ 
ing  the  low-held  approximation,  which  shows  increasingly 
severe  inaccuracies  as  i/io  exceeds  1,  especially  when  the 
symmetry  factor  deviates  from  0.5  in  either  direction.  Fig.  5 
shows  significant  errors  using  the  high-held  approximation 
for  i/i'o<4  and  £ < 0.5.  The  range  of  applicability  increases 
for  /I  >0.5.  Fig.  6  indicates  small  error  over  the  entire  op¬ 
erating  regime  using  the  hyperbolic  sine  approximation  for 
£<0.7. 


i/i„ 


4.  Combining  anode  and  cathode  activation  losses 

Another  questionable  practice  in  some  SOFC  models  is  to 
represent  the  combined  anode  and  cathode  activation  losses 
with  a  single-term  equation.  For  example,  starting  with  the 
hyperbolic  sine  approximation,  the  total  cell  activation  loss 
is: 


A  Vact, total  =  AVactiA+  A  Vact, c 


There  are  two  ways  to  reduce  this  two-term  equation  down 
to  the  following  single-term  equation  and  retain  the  physical 
meanings  of  a  and  ;q: 


Fig.  5.  Deviation  from  B-V  equation  using  high-field  approximation. 


(15) 
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Method  1:  Assume  identical  activation  loss  behavior  for 
the  anode  and  cathode  (i.e.,  «i;a  =«i,c  and  to, a  =  io,c,  which 
yields  aeff  =  o!i,a/2  and  i'o,eff=io,A)-  This  method  is  ex¬ 
tremely  questionable,  since  the  activation  loss  is  dependent 
on  electrode  materials  and  local  reactants  and  products,  which 
are  completely  different  for  the  anode  and  cathode  of  a 
SOFC. 

Method  2:  Assume  one  of  the  equilibrium  exchange  cur¬ 
rent  densities  (i.e.,  anode  or  cathode)  is  sufficiently  larger  than 
the  other,  thus  allowing  the  corresponding  activation  loss  to 
be  neglected  (e.g.,  i(i\  >  1  OOio.c  is  required  to  neglect  the  an¬ 
ode  activation  loss  and  keep  the  resulting  error  below  2%).  If 
this  is  not  the  case  and  experimental  data  is  used  to  develop 
curve-fits  for  values  of  aeff  and  *o,eff,  then  the  model’s  pre¬ 
dictive  capabilities  should  be  considered  suspect  outside  the 
operating  conditions  of  the  test  unit. 


5.  Conclusions 

For  the  relatively  low  temperature  SOFCs  being  devel¬ 
oped  today,  it  is  no  longer  reasonable  to  neglect  the  ac¬ 
tivation  losses  when  building  models  to  predict  the  cell 
voltage-current  characteristics.  The  Butler- Volmer  equation, 
or  an  appropriate  approximation,  should  be  implemented  to 
accurately  predict  these  losses.  If  an  explicit  equation  is  de¬ 
sired  for  functional  reasons  (e.g.,  quick  parametric  studies 
or  severe  computational  speed  requirements)  and  the  charge 
transfer  reaction  is  assumed  to  be  a  one-step,  single-electron 
transfer  process,  then  the  hyperbolic  sine  approximation  (Eq. 
(13))  is  recommended  for  its  superior  accuracy  over  the  ma¬ 
jority  of  operating  current  densities  and  symmetry  factors 
of  interest.  Otherwise,  the  implicit  Butler- Volmer  equation 
is  easily  solved  using  any  one  of  a  variety  of  standard  nu¬ 
merical  technique  (e.g.,  the  basic  secant  method  typically 
converges  within  five  iterations).  High-held  and  low-held  ap¬ 
proximations,  such  as  Eqs.  (11)  and  (12),  should  be  avoided 
unless  the  model’s  range  of  i/io  values  is  clearly  established  to 
be  within  the  approximation’s  range  of  applicability.  Finally, 
if  a  single-term  equation  is  used  to  combine  the  anode  and 
cathode  activation  losses,  then  the  applicability  of  this  model 
should  be  scrutinized  outside  of  the  range  of  the  experimental 
data  used  for  validation. 
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